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1 



1 "Borehole Surveying" 
2 

3 This invention relates to a method and apparatus for 

4 use in surveying of boreholes. 
5 

6 It is known in directional drilling, for example, to 

7 detect the orientation of a drillstring adjacent to 

8 the bit by means of a sensor package for determining 

9 the local gravitational [GX,GY,GZ] and magnetic 

10 [BX,BY,BZ] field components along mutually 

11 orthogonal axes, and to derive from these the local 

12 azimuth (AZ) and inclination (INC) of the 

13 drillstring. Conventionally, the measurements are 

14 made by providing within the instrument package 

1 5 three mutually perpendicular acceleromet ers and 

16 three mutually perpendicular magnetic f luxgates . 
1 7 

18 The present invention is concerned with an 

1 9 arrangement which requires only two measurement 

20 devices, namely a single accelerome ter and a single 

21 magnetic fluxgate or a single acceleromet er and a 

22 single rate gyro, the latter being preferred for 

23 situations in which magnetic interference is likely 

24 to be encountered. 
25 

26 Accordingly, the present invention provides a method 

27 of surveying boreholes, comprising: 



2 

1 providing an instrument package in the leading 

2 end of a drillstring, the instrument package 

3 comprising first and second single-axis sensors 

4 mounted for rotation with the drillstring about the 

5 rotational axis of the drillstring, the first sensor 

6 being an accelerometer and the second sensor being a 

7 magnetic fluxgate or a rate gyro; 

8 rotating the drillstring; 

9 deriving from the first sensor the inclination 
10 angle of the drillstring at the instrument package; 
1 1 and 

12 deriving from the second sensor the azimuth 

13 angle of the drillstring at the instrument package. 
1 4 

15 Each of the sensors will typically be positioned in 

16 one of two configurations. In the first 

17 configuration, the sensor is radially spaced from 

18 the borehole axis and has its sensing axis in a 

19 plane containing the borehole axis and an axis 

20 perpendicular thereto. In the second configuration, 

21 the sensor is radially spaced from the borehole axis 

22 and has its sensing axis in a plane parallel with 

23 the borehole axis. 
24 

25 Preferably, the drilling control rotation angle is 

26 also obtained from the sensor outputs. 
27 

28 Preferably, the sensor outputs are integrated over 

29 the four quadrants of rotation and the desired 

30 output angle is derived from the integrated output. 

31 The instrument package suitably includes rotation 

32 angle reference means for use in the integration. 



3 



1 

2 Additional information may be derived, such as the 

3 local gravitational and magnetic field vectors. 
4 

5 From another aspect, the invention provides 

6 apparatus for use in surveying boreholes, the 

7 apparatus comprising an instrument package adapted 

8 to be included in the leading end of a drillstring, 

9 the instrument package comprising first and second 

10 single-axis sensors mounted for rotation with the 

1 1 drillstring about the rotational axis of the 

12 drillstring, the first sensor being an accelerometer 

1 3 and the second sensor being a magnetic f luxgate or a 

14 rate gyro; and computing means for deriving from the 

1 5 first sensor while the drillstring is rotating the 

1 6 inclination angle of the drillstring at the 

17 instrument package, and for deriving from the second 

1 8 sensor while the drillstring is rotating the azimuth 

19 angle of the drillstring at the instrument package. 
20 

21 The computing means preferably operates to integrate 

22 the sensor outputs over the four quadrants of 

23 rotation and to derive the desired output angle 

24 from the integrated output. 
25 

26 The apparatus may further include rotation angle 

27 reference means for use in the integration. 
28 

29 Examples of the present invention will now be 

30 described, by way of illustration only, with 

31 reference to the drawings, in which: 
32 



4 

1 Fig. 1 illustrates, in general terms, the 

2 operation of a single axis sensor in a drillstring 

3 for sensing any given vector V; 

4 Fig. 2 is a block diagram of one circuit which 

5 may be used to identify rotation quadrant; 

6 Fig. 3 illustrates the operation where the 

7 sensor is an acceleromet er ; 

8 Fig. 4 illustrates the operation where the 

9 sensor is a fluxgate; 

10 Fig. 5 illustrates the derivation of azimuth 

1 1 angle ; and 

12 Fig. 6 illustrates the operation where the 

1 3 sensor is a rate gyro . 
1 4 

1 5 

1 6 Single -axis sensor 

18 The operation of a single-axis sensor in a drill 

19 string will first be described in general terms. 

20 The application of this to specific sensors is 

21 discussed below. 
22 

23 Referring to Fig. 1, a single-axis sensor 10 is 

24 mounted on a drill string (not shown) . The sensor 

25 10 senses a fixed vector {V} and is mounted in one 

26 of two configurations. 
27 

28 In the first configuration, the sensor 10 lies in a 

29 plane containing the rotation axis (OZ) of the drill 

30 string and axis (OX) perpendicular to (OZ) . Axis 

31 (OY) makes up the conventional orthogonal set of 

32 axes [OX,OY,OZ]. The sensor 10 is mounted at a 



5 

1 distance r from the (OZ) axis and the angle between 

2 the sensing axis (OS) and the rotational axis (OZ) 

3 ism. 

4 

5 In the second configuration, the sensor 10 is 

6 mounted in a plane which is parallel to the borehole 

7 axis (OZ) and with its sensing axis perpendicular to 

8 the axis (OY) and making angle m with the direction 

9 of the borehole axis (OZ) . 

I 0 

II If the rate of rotation about the (OZ) axis is w and 

12 the components of {V} are {VOZ} along the (OZ) axis 

13 direction and {VOXY} in the (OXY) plane, then if the 

14 output from the sensor 10 for both configuration 1 

1 5 and configuration 2 of Figure 1 is of the form 
1 6 

17 V(t) = VOZ.cos(m) + VOXY . sin (m) . cos (w . t ) + c 

1 8 

19 where time t = 0 when the axis (OX) is coincident 

20 with the direction of {VOXY} and c is constant for 

21 any fixed rotation rate w. 
22 

23 Thus, the sensor output at time t can be written: 
24 

25 V(t) = Kl.cos(w.t) + K2 (i) 

26 

27 where K1 = VOXY . sin (m) and K2 = VOZ. cos (m) + c are 

28 constant if the vector amplitudes VOZ and VOXY are 

29 constant. 
30 
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1 Sensor output integration 
2 

3 The integration of V(t) from any initial time ti to 

4 ti + T/4, where T = 2.3t/w, the time for one 

5 revolution about (OZ) , is 
6 

H+T/4 ti+T/4 

7 Q= \ K I. cos( wx). dt + \K2.dt 
8 

9 Thus, 

10 ti + T/4 

11 Q= [ (K1/w) .sin(w.t) ] +K2.T/4 

1 2 ti 

1 3 

1 4 or 
1 5 

16 Q = (K1/w) . [sin(w.ti + w.T/4) - sin(w.ti)] + L 

1 7 

1 8 or 

19 Q = (K1/w) . [sin(w.ti + 3t/2) - sin(w.ti)] + L 

20 or 

21 Q = (K1/w) . [cos (w.ti) - sin(w.ti)] + L (ii) 

22 where L is a constant = K2.T/4. 
23 

24 Using equation (ii), the integration of V(t) from an 

25 arbitrary time tO to time tO+T/4 yields 
26 

27 Q1 = (K1/w) . [cos(w.to) - sin(w.to)] + L (iii) 

28 

29 Using equation (ii), the integration of V(t) from 

30 time tO+T/4 to time tO+T/2 yields 
31 



7 

1 Q2 = (K1/w) . [cos(w.tO +W.T/4) - sin(w.tO + w.T/4)]+L 

2 or 

3 Q2 = (K1 /w) . [cos (w. tO + Jt/2) - sin(w.tO + jr/2)]+L 

4 or 

5 Q2 = (K1/w) . [-sin(w.tO) - cos(w.tO)] + L ...(iv) 

6 

7 Using equation (ii), the integration of V(t) from 

8 time tO+T/2 to tO+3T/4 yields 
9 

10 Q3 = (K1/w) . [cos (w. tO+w.T/2) - sin ( w . tO+w . T/2 ) ] +L 

1 1 or 

12 Q3 = (K1 /w) . [cos (w. tO+jt) - sin (w. tO+rt ) ] + L 

1 3 or 

14 Q3 = (K1/w) . [-cos(w.tO) + sin(w.tO)] + L (v) 

1 5 

16 Using equation (ii), the integration of V(t) from 

17 time tO+3T/4 to time tO+T yields 
1 8 

19 Q4 = (K1 /w) - [cos (w.tO+w.3T/4) - sin (w. tO+w. 3T/4) ] +L 

20 or 

21 Q4 = (K1 /w) . [cos (w. t0+3jt/2) - sin(w. t0+3jt/2) ] +L 

22 or 

23 Q4 = K1/w) . [sin(w.tO) + cos(w.tO)} + L (vi) 

24 

25 Writing K = K1 /w and a = w.tO, then equations (iii) 

26 through (vi) yield for the four successive 

27 integrations of V(t) 
28 

29 Q1 = -K.sina+ K.cosa + L (vii) 

30 Q2 = -K.sina - K.cosa + L (viii) 

31 Q3 = K.sina - K.cosa + L (ix) 



8 

1 Q4 = K. since + K.cosa + L (x) 

2 

3 Integration control 
4 

5 In order to control the sensor output integration, 

6 as just described, over four successive quarter 

7 periods of the drill string rotation, a train of n 

8 (with n any multiple of 4) equally spaced pulses per 

9 revolution must be generated. If one pulse P 0 of 

10 this pulse train is arbitrarily chosen at some time 

11 tO, the repeated pulses P n/4/ P n /2 and P 3n /4 define 

12 times tO+T/4, tO+T/2 and tO+3T/4 respectively where 
1 3 the period of rotation T = 2jt/w and w is the angular 
14 velocity of rotation. 

1 5 

16 A suitable means for generating an appropriate 

1 7 control pulse train is described in US-A1 - 

18 20020078745, which is hereby incorporated by 

1 9 reference . 
20 

21 In an alternative form of integration control, the 

22 sensor output waveform itself can be used with 

23 appropriate circuitry for defining the integration 

24 quadrant periods. In particular, the relatively low 

25 noise magnetic fluxgate output is well suited to act 

26 as input to a phase-locked-loop arrangement. Fig. 2 

27 shows such an arrangement, successive output pulses 

28 defining the integration quadrants. 
29 

3 0 Rotation angle 
31 
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1 Equations (vii) through (x) can be solved to yield 

2 angle a; there is a degree of redundancy in the 

3 possible solutions but, for example, 
4 

5 Q1 - Q2 = 2K.cosa 

6 and 

7 Q3 - Q2 = 2K. since 

8 or 

9 sina/cosa = (Q3-Q2 ) / (Q1 -Q2 ) (xi) 

I 0 

II Since a= w.tO, the angle S(tO) between the axis 

12 (OX) and the direction of {VOXY} at time tO can be 

13 determined from equation (xi), and the angle between 

14 (OX) and {VOXY} at any time tm measured from the 

1 5 arbitrary starting time tO is then 
1 6 

17 S(tm) = a + w.tm = S(tO) + 2jt.tm/T ....(xii) 

1 8 

1 9 Magnitudes of vectors {VOXY} and {VOZ} 

20 

21 Equations (vii) through (x) can be solved to yield 

22 the constant L: 
23 

24 L = (Q1 + Q2 + Q3 + Q4)/4 (xiii) 

25 

26 and the constant K can be determined from: 
27 

28 (K) 2 = [(Q1-L) 2 + (Q2-L) 2 ]/2 

29 = [(Q3-L) 2 +(Q4-L) 2 ]/2 . . . (xiv) 

30 

31 The magnitude of vector {VOZ} can be determined as 



1 0 



1 

2 VOZ = (K2-c) /cos (m) = (4.L/T - c)/cos(m) . . . . (xv) 

3 provided that constant c is known. 
4 

5 The magnitude of vector {VOXY} can be determined as 
6 

7 VOXY = K1/sin(m) = (K.w) /sin(m) (xvi) 

8 

9 Inclination angle 

I 0 

II The inclination angle (INC) can be derived from the 
1 2 gravity vector {G} with the aid of a rotating 

1 3 accelerometer . 
14 

15 Referring to Fig. 3, where (INC) is the angle 

16 between the tool axis (OZ) and the gravity vector 

17 {G}, 
1 8 

19 GOZ = G.cos(INC) (xvii) 

20 and 

21 GOXY = -Gsin(INC) (xviii) 

22 

23 The accelerometer output can be written as 
24 

25 VG(t) = GOZ.cos(m) + GOXY . sin (m) . cos (wt ) 

26 + CP.sin(m) + D.sin(m) (xix) 

27 

28 where CP is a centripetal acceleration term and D is 

29 a sensor datum term. The centripetal acceleration 

30 term CP is zero for configuration 2 and makes this 

31 the preferred configuration for mounting of the 

32 accelerometer. 



1 1 

1 

2 Since CP is proportional to w 2 /r and is constant for 

3 constant w, then clearly VG(t) is of the form 
4 

5 VG(t) = Kl.cos(w.t) + K2(w) 

6 (or Kl.cos(w.t) + K2 for configuration 2) (xx) 

7 

8 where K1 and K2(w) are constants at constant angular 

9 velocity w in the case of configuration 1 and always 

10 constant in the case of configuration 2. the 

11 constants K1 and K2 (w) can be determined from the 

12 accelerometer output integrations as described above 

13 together with the angle (Highside Angle HS = w.t) 

14 between the axis (OX) and the direction of {GOXY} . 
1 5 

16 K1 = GOXY.sin(m) (xxi) 

1 7 and 

18 K2(w) = GOZ.cos(m) + D.sin(m) (xxii) 

19 with 

20 C(w) = CP.sin(m) + D.sin(m) (xxiii) 

21 constant at constant angular velocity w (or for 

22 configuration 2 at all w) . 
23 

24 A calibration procedure can be carried out to 

25 determine the values of C(w) for angular velocity 

26 values w (constant in the case of configuration 2) 

27 by calculating values of K2 (w) with the rotation 

28 axis (OZ) horizontal when C(w) = K2(w). 
29 

30 Thus, for any drilling situation with known angular 

31 velocity w, the vector components of the local 

32 gravity vector {G} can be determined as 



1 2 



1 

2 GOXY = K1/sin(m) (xxiv) 

3 and 

4 GOZ = (K2(w) - C(w))/cos(m) (xxv) 

5 



6 The inclination angle (INC) can then be determined 

7 from 
8 

9 sin (INC) /cos (INC) = -GOXY/GOZ (xxvi) 

1 0 

1 1 Azimuth angle 
1 2 

13 When using a rotating fluxgate, the azimuth angle 

14 (AZ) can be determined from a consideration of the 

15 magnetic vector {B} . What follows is applicable to 

16 both configuration 1 and configuration 2. 
1 7 

18 With reference to Fig. 4, it can be shown that 

1 9 

20 BOZ = BV. cos (INC) 

21 + BN.cos (AZ) . sin(INC) . (xxvii) 

22 

2 3 and 
24 

25 BOXY = (BN.cos (AZ) .cos (INC) -BV.sin(INC) ) .cos (HS-MS) 

26 + BN.sin(AZ) .sin(HS-MS) (xxviii) 

27 

28 or, with HS-MS = d a constant, 
29 

30 BOXY = (BN.cos(AZ) . cos ( INC) -BV . sin ( INC) ) .cos(d) 

31 +BN.sin(AZ) .sin(d) (xxix) 

32 



1 3 

With D the fluxgate datum, the fluxgate output can 
be written 



VB(t) = BOZ.cos(m) + BOXY.sin(m) .cos(w.t) 
+ D.sin(m) 



(xxx) 



or 



VB(t) = Kl.cos(w.t) + K2 



(xxxi) 



where 



K1 = BOXY.sin(m) 



and 



K2 = BOZ.cos(m) + D.sin(m) 



= BOZ.cos(m) + C 



(xxxii) 



are constants which can be determined from the 
fluxgate output integrations as described above 
together with the angle (Magnetic Steering Angle = 
MS = w.t) between the axis (OX) and the direction of 
{BOXY} . 

A calibration procedure can be carried out to 
determine the value of the constant C by calculating 
the value of K2 while rotating about the direction 
of the axis (OZ) along which BOZ = 0 when K2 = C. 

Thus, for any drilling situation the vector 

components of the local magnetic field {B} can be 
determined as 

BOXY = K1/sin(m) (xxxiii) 



and 



BOZ = (K2-C) /cos (m) 



(xxxiv) 



1 4 

1 With reference to Fig. 5, the horizontal component 

2 {BN} of the local magnetic field vector {B} can be 

3 represented by horizontal components {B1 } and {B2} 

4 where 
5 

6 B1 = BOXY-cos (d) .cos (INC) 

7 + BOZ.sin(INC) (xxxv) 

8 and 

9 B2 = BOXY.sin(d) (xxxvi) 

I 0 

I I The Azimuth Angle (AZ) can then be determined from 
1 2 

13 sin(AZ) /cos (AZ) = -B2/B1 (xxxvii) 

1 4 

15 Also, the horizontal component of the local magnetic 

1 6 field can be determined from 
1 7 

18 BN = (B1 2 + B2 2 ) 3/2 (xxxviii) 

1 9 

20 and the vertical component of the local magnetic 

21 field can be determined from 
22 

23 BV = BOZ.cos(INC) 

24 - BOXY. cos (d) .sin (INC) (xxxix) 

25 

2 6 Earth's rotation vector 
27 

28 Where it is not practicable to use a magnetic 

29 fluxgate, this may be replaced by a rate gyro as 

30 sensor. 
31 



1 5 

1 With reference to Fig. 6, if the geographic latitude 

2 at the drilling location is (LAT) then the vertical 

3 component of the earth's Rotation Vector {RE} is 
4 

5 RV = -RE . sin (LAT) (xl) 

6 and the horizontal component is 

7 RN = RE.cos(LAT) (xli) 

8 

9 The magnitude of the cross-axis rate vector {ROXY} 

1 0 can be shown to be 
1 1 

12 ROXY = (RN.cos (GAZ) .cos (INC) -RV.sin(INC) ) .cos (d) 

13 + RN.sin(GAZ) sin(d) (xlii) 

14 

15 where (GAZ) is the gyro azimuth angle and 

16 d = HS - GS is constant. 
1 7 

18 Since RN, RV, d and INC are known and ROXY can be 

19 derived as discussed below, (GAZ) can be determined. 
20 

21 With the particular configuration where the rate 

22 gyro sensing axis is perpendicular to the drill 

23 string rotation axis (OZ) , the rate gyro output can 

24 be written 
25 

26 VG(t) = ROXY. cos ( w. t ) + D (xliii) 

27 

28 where D is the rate gyro datum, or 
29 

30 VG(t) = Kl.cos(w.t) + K2 (xliv) 

31 



1 6 

1 where the constant K1 = ROXY can be determined from 

2 the rate gyro output integrations as described above 

3 together with the Gyro Steering Angle GS = w.t 

4 between (OX) and the direction of {ROXY}. 
5 

6 The variation in the Rate Gyro Datum makes it 

7 difficult to achieve satisfactory datum calibration 

8 in all circumstances. It is unlikely that Gyro 

9 Azimuth measurements should be attempted at high 

10 inclination angles. The use of the rate gyro is 

1 1 most likely with near-vertical boreholes in 

1 2 locations where magnetic azimuth measurements are 

13 unreliable (such as close to rigs) and the Gyro 

14 Azimuth GAZ is approximately equal to the angle d. 
1 5 

16 The present invention thus makes possible the 

17 measurement of a number of borehole-related 

18 parameters during rotation of a drillstring and 

19 using a reduced number of sensors. Modifications 

20 may be made to the foregoing embodiments within the 

21 scope of the present invention. 



